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Abstract 

It is well known that the effect of gamma-ray absorption on extragalactic back¬ 
ground light (EBL) is weakly expressed in the spectra of some blazars. It is shown 
that a secondary component generated by electromagnetic cascades might consid¬ 
erably decrease the statistical significance of this anomaly. Observational results 
indicate the existence of the cascade component in the spectra of extragalactic 
ganrma-ray sources, thus supporting the proposed model. 

Keywords: blazars: Mkn 421, Mkn 501, H 1426+428, 1ES 1101-232, 1ES 0347-121, 
1ES 0414+009; extragalactic background light (EBL); 77 absorption; EM cascades; ex¬ 
tragalactic magnetic field (EGMF). 

1 Introduction 

For a long time it was expected that gamma-ray spectrum of a distant source in the 
very high energy (VHE, E >100 GeV) range would have a steepening resulting from the 
77 —> e + e~ absorption n-H on the extragalactic background light (EBL) photons at the 
energy region where the optical depth of this process t 77 = 1. However, for some blazars 
- gamma-ray loud active galactic nuclei — the ratio of observed to expected intensity 
K after accounting for the EBL absorption effect rises with r 77 ; K is larger for r 77 > 2 
than for 1 < r 77 < 2 [3]. The statistical significance of this anomaly is 4.2 a [3]. 

An effect of such kind is under discussion for about 15 years 0 0 ; during this period 
of time two popular models of the anomaly were developed. The first model assumes 
that blazars are sources of hadrons that produce photons on the way from the source 
to the observer i-0- I 11 this case a good collimation of the hadronic beam is required 
(the same as for the gamma-ray beam, i.e. 9j et ~ 1 °); this condition is hard to meet, as 
after the exit from the blazar’s jet a region of chaotic magnetic field is likely to exist, and 
this magnetic field is able to scatter the accelerated hadrons. The second model (e.g., 
0) postulates oscillations of photons to light exotic particles in magnetic field (and vice 
versa); the efficiency of the reverse conversion to photons is required to be quite high, 
independently of the direction to the source, and, therefore, of the Galactic magnetic held 
parameters in this direction. 

Another possible mechanism that could explain deviation of K in the optically thick 
regime from the expected value is production of secondary particles in electromagnetic 
cascades a. i-m- The present work is devoted to examination of the cascade compo¬ 
nent’ influence to the statistical significance of the discussed anomaly. 
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2 Observational data and analysis methods 

The present work is based on the sample of objects that were observed by Cherenkov 
telescopes listed in [3], The following spectra were selected, that contain: 1) at least 6 
bins 2) at least one bin in the r 77 > 2 region, in total 6 objects: Mkn 421 [11], Mkn 501 
g], H 1426+428 pj], 1ES 1101-232 g3], 1ES 0347-121 [H], 1ES 0414+009 [[5] (redshifts 
0 from 0.031 to 0.287). 

The full MC calculation of ~ 10 6 cascade spectra from primary photons in the energy 
range E 0 = 100 GeV-100 TeV was performed with the ELMAG 2.02 code [T 6 ], assuming 
the EBL model of m- The shape of the primary spectrum was set to 
J = C ■ E 0 7 • exp(—Eo/E c )i where C is some constant, 7 — the primary spectrum index, 
E c — the energy corresponding to the end region of the spectrum. The ID assumption 
used in the calculations is well justified, as in the energy range under investigation the 
angular radius of a typical cascade 9 C « 0 Jet . 

Various hypotheses were tested for the described sample of objects; for each hypoth¬ 
esis the corresponding p-value was calculated and converted to the equivalent statistical 
significance (Z-value, see [18], Chapter 2). For the significance calculation the values of 
the observed and model intensity in certain ranges of the r 77 value were utilized ((l 8 j. 
Chapter 5.1). It was assumed that the distributions of the measured values in both di¬ 
rections (increasing and decreasing one) are Gaussian, but with the different parameters; 
a more comprehensive analysis calls for some more general, asymmetric distribution (e.g. 

ra-ra). 

3 The anomaly without the cascade component 

If cascade electrons are deflected by the extragalactic magnetic field (EGMF) to the angles 
9m » 9j et , the cascade component in the observed spectrum is practically absent. An 
example of the spectral energy distribution (SED, E 2 dN/dE ) for one of the analysed 
spectra in this case is shown in Fig. [H The model of the observed spectrum (green line) 
doesn’t contradict the observations at the r 77 < 2 region, but for r 77 > 2 the model 
intensity is somewhat lower than needed to desribe the data. 

The hypothesis that without the cascade component the model is still able to describe 
the data in the r 77 > 2 region for the sample of 6 objects was rejected at the statistical 
significance level Z\ = 3.1cr. A similar result was obtained in [3] with a different method: 
Z = 3.8 a. I 11 general, these results are in good coincidence, because the systematic 
uncertainty of the spectral shape that was accounted for in the present work tends to 
decrease the significance. 

4 Influence of the cascade component to the statisti¬ 
cal significance of the anomaly 

For the EGMF strength B <1CT 16 G and correlation length L c = 1 Mpc, cascade photons 
can significantly contribute to observed spectrum at E > 100 GeV (see Fig. [2j). While the 
77 —» e + e _ process due to its threshold nature occurs only at EBL photons for F 7 <100 
TeV, the main target for electrons in the Inverse Compton (IC) interaction acts is dense 
cosmic microwave background (CMB) [9], [16], whose photons have comparatively low 
energy. Therefore, the IC process, as a rule, occurs in the Thomson regime, thus the 
photons of the next generation carry only a small fraction of the primary electron energy. 
This peculiarity of the cascade process is clearly seen in Fig. [2} the secondary component 
from cascade photons (blue curve) gives the main contribution at comparatively low 
energy, in the optically thin regime. 

In the framework of the presented model, the contribution of the cascade component 
at t 77 > 1 cannot account for the anomaly in the distant (z > 0 . 2 ) blazar spectra at 
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Figure 1: Illustration of a spectral energy distribution analysis for 1ES 1101-232 (z= 
0.186) [IS] without the cascade component. Red circles — the measured spectrum, red 

dashed line — the total uncertainty of the measured spectrum yja% at + ^1 yst . with account 

of statistical cr staL and systematic a syst . uncertainties, black line — the primary model 
spectrum, green line — the model of the observed spectrum, blue vertical line — the 
energy, at which t 77 = 1, blue vertical line — the same, but for t 77 = 2. 


E > 1 TeV. However, the secondary photons do influence the interpretation of the data 
indirectly: for better agreement of the model with observations, the primary spectrum in 
Fig. [2] is chosen harder, than in Fig. [TJ this circumstance, in its turn, leads to a significant 
increase of the model intensity at the r 77 > 2 region. An analysis similar to the one 
conducted in sec. |3l but with account of the cascade component, yields Za= 0.36 cr (the 
anomaly is practically absent). 


5 Discussion 

The analysis of the cascade component influence on the gamma-ray spectrum in the 
present work was carried out assuming B <<10” 16 G. The probable influence of the 
EGMF on the cascade component spectrum is schematically shown by dashed red curve 
in Fig. [2j A similar spectral shape was already observed during a flare of Mkn 501 [21] , that 
was simultaneously observed by the orbital Fermi LAT instrument [22] and the VERITAS 
Cherenkov telescope [23]. Additionally, in [21] a constraint on the impulse front rise time 
Tp <10 days (90 % C.L.) at E~100 GeV was obtained, as well as an estimate of B= 
10 -16 -10 -17 G for L c = 1 Mpc. 

If B and L c are nearly the same for all 6 objects of the considered sample, the char¬ 
acteristic variability time of the cascade component for them could be estimated as [3] 
{zf z(Mkn501)) ■ T F = 10-80 days, depending on z. For the two nearest objects of the 
sample the contribution of the cascade component is small; 4 blazars with z > 0.1 do not 
show variability with period less than several months. This doesn’t allow to set additional 
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Figure 2: Analysis of the same spectrum as in Fig. Q] but with account of the cascade 
component. The notations are the same as in Fig. |TJ except: 1. blue curve that denotes 
the model spectrum with account of the cascade component 2. red dashed curve that 
qualitatively illustrates a possible influence of EGMF on the cascade component spectrum. 


constraints on the cascade component’ contribution using temporal information. 

Existing bounds on B (e.g., m, where the angular distribution of observed photons, 
that depends on the cascade component contribution, is analysed) doesn’t contradict the 
assumptions of the considered model. Future instruments with either high sensitivity (e.g. 
CTA [25]) or extremely good angular resolution (e.g. emulsion gamma-ray telescope paj 
or the GAMMA-400 apparatus m) would possibly allow to measure the EGMF strength. 
Finally, the recent indication that the directions to the objects observed by Cherenkov 
telescopes and the Fermi LAT instrument are predominantly point to large-scale voids 
[28], and, therefore, to low B regions, is a direct argument in favor of existence of the 
cascade component in blazar spectra. 

6 Conclusions 

In the present work, for the first time, a possibility to explain the anomaly in blazar 
spectra at VHE region in the framework of cascade model is considered in detail (i.e. 
quantitatively). For the first time, it is shown that the influence of secondary photons 
that might decrease the statistical significance of the anomaly may have an indirect nature: 
while these photons contribute mainly in the optically thin regime, the cascade component 
allows the harder primary spectrum to fit the data, which leads to better agreement of 
the model with observations in the optically thick regime. 
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The cascade model doesn’t require a single additional assumption and doesn’t contra¬ 
dict to observations of Cherenkov and orbital gamma-ray telescopes; at the present time 
this model is the simplest explanation of the anomaly for some blazars. 

Additional note. An evidence for broadened angular distribution at ~1 GeV for a sam¬ 
ple of low redshift blazars was recently found in I2®j using Fermi LAT data and interpreted 
as a proof of pair halos existence; an estimate of EGMF strength of 1(F 15 ~1(F 17 G was 
obtained. The latter value doesn’t contradict the model presented here. 
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